Optimal discrimination of mixed quantum states involving inconclusive results 
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We propose a generalized discrimination scheme for mixed quantum states. In the present scenario 
we allow for certain fixed fraction of inconclusive results and we maximize the success rate of the 
quantum-state discrimination. This protocol interpolates between the Ivanovic-Dieks-Peres scheme 
and the Helstrom one. We formulate the extremal equations for the optimal positive operator valued 
measure describing the discrimination device and establish a criterion for its optimality. We also 
devise a numerical method for efficient solving of these extremal equations. 
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I. INTRODUCTION 



The non-orthogonality of quantum states is one of the 
fundamental features of quantum mechanics. It means 
that two different quantum states cannot be perfectly 
discriminated in general. This fact imposes limits on 
the processing of quantum information. For instance, 
an unknown quantum state cannot be copied ^ |^. On 
the other hand, the non-orthogonality of the quantum 
states can be advantageous. In particular, most of the 
quantum cryptographic schemes rely on it H, IH || . The 
principle of operation of these schemes is to employ sev- 
eral non-orthogonal states to encode the bits of the se- 
cret key that is sent from transmitter (Alice) to receiver 
(Bob). Since the eavesdropper (Eve) cannot discriminate 
all these states perfectly, any attempt to gain information 
on the key will increase the noise of the transmission from 
Alice to Bob and thus Eve will be unavoidably revealed. 

Since one cannot discriminate non-orthogonal states 
perfectly, a natural question arises what is the opti- 
mal approximate discrimination scheme. Two different 
approaches to the problem have been proposed. The 
first scenario, the ambiguous quantum state discrimina- 
tion, has been analyzed by Holevo, Helstrom and others 
|, 0, |, |, |10|, |l^, |l|, 0. They considered discrimina- 
tion of mixed quantum states pj that are generated 
by Alice with the a-priori probabilities pj , 
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pj 



Pj > 0. 



(1) 



Bob performs a generalized quantum measurement de- 
scribed by the A^-component positive operator valued 
measure (POVM) {Iljlj^j^. If the outcome Hj is de- 
tected, then Bob concludes that Alice sent the state pj 
to him. The optimal POVM is defined as the POVM that 
maximizes the average success rate of the Bob's guesses, 



N 



^s-EP^Tr[n,p,]. 



(2) 



An interesting alternative approach has been suggested 
by Ivanovic, Dieks, and Peres (IDP) for the discrimina- 
tion of two pure states ^ |l^, |l^ and extended to 



A^ linearly independent states by Chefles and Barnett 
jTs], |l9[ ^ 0. These states can be discriminated un- 
ambiguously, provided that we allow for some fraction 
of the inconclusive results Pi. Recently, several other 
problems closely related to quantum-state discrimination 
have been considered such as the discrimination of sets 
of quantum states ||2^, ^ and the quantum-state com- 
parison p5| . 

Besides being a theoretically interesting problem, the 
quantum-state discrimination has also found practical 
applications in quantum information processing. In par- 
ticular, the optimal unambiguous quantum-state discrim- 
ination represents a simple and experimentally feasible 
attack on the quantum key distribution protocols |2^, ^ . 
Recently, some optimal discrimination POVMs have been 
realized experimentally and the Helstrom as well as the 
IDP bounds have been attained [|8[|9[ |o[ llj. 

The two above discussed scenarios can be considered as 
limiting cases of a more general scheme that involves cer- 
tain fraction of inconclusive results Pi for which we max- 
imize the success rate. The main feature of this scheme 
is that if we allow for inconclusive results then we can 
improve the relative (or re-normalized) success rate 
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1 -Pt 



(3) 



In other words, with probability Pj Bob fails completely 
and he cannot say at all which state was sent to him. 
However, in the rest of the cases when he succeeds he can 
correctly guess the state with higher probability than if 
he would not allow for the inconclusive results. For pure 
linearly independent states this generalized scenario was 
discussed in two recent papers ^ . 

Here we extend the analysis to general mixed states. 
At a first sight, such an extension may seem to be a bit 
problematic, because it is known that one cannot unam- 
biguously discriminate mixed states (the reason is that 
the IDP scheme does not work for linearly dependent 
states). Nevertheless, we shall show that the extension is 
perfectly meaningful also for mixed states. Although we 
cannot reach the limit of perfect unambiguous discrimi- 
nation, we can improve the Prs and we shall derive an 
upper bound on Prs that can be achieved for a given set 
of mixed states. 
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II. EXTREMAL EQUATIONS FOR OPTIMAL 
POVM 



Let us begin with the formal definition of the problem. 
We assume that the quantum state sent to Bob is drawn 
from the set of N mixed states {pjjjLi with the a-priori 
probabilities pj. Bob's measurement on the state may 
yield + 1 different results and is formally described by 
the POVM whose iV + 1 components satisfy 



N 

En. 



= 1, 



(4) 



where 1 is the identity operator. The outcome Hq indi- 
cates failure and the probability of inconclusive results is 
thus given by 
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(5) 



For a certain fixed value of Pi we want to maximize the 
relative success rate ^ which is equivalent to the maxi- 
mization of the success rate . To account for the linear 
constraints (Q) and (|^) we introduce Lagrange multipli- 
ers A and a where A is Hermitian operator and a is a real 
number. Taking everything together we should maximize 
the constrained success rate functional 
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3 = 1 



N 

E 

j=0 
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Tr[An, 



j=i 



(6) 

We now derive the extremal equations that must be sat- 
isfied by the optimal POVM. We expand the POVM 
elements in terms of their eigenstates and eigenvalues, 
= J2k^3k\fjk){V3k\ and vary (§) with respect to 
{^jk\- After some algebraic manipulations, we arrive at 
the extremal equations, 

(A-pjPj)n, =0, j = l,...,iV, (7) 
(A-acr)no = 0, (8) 

where the operator a introduced for the sake of notational 
simplicity reads 
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(9) 



From the constraint Tr[ano] — Pi we can express a in 
terms of A, 

a = PfiTr[Ano]. (10) 

Furthermore, if we sum all Eqs. @) and also Eq. (|) and 
use the resolution of the identity we obtain formula 
for A, 
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A = ^pj-pjEj + flfTlIo. 



(11) 



If we combine Eqs. and then we can express 

a and A solely in terms of pj , pj , and Hj . This may be 
important, for example, if we guess the optimal POVM 
and want to determine the corresponding Lagrange mul- 
tipliers. The extremal Eqs. (j^) and (||) constitute a gen- 
eralization of the extremal equations for optimal POVM 
for ambiguous quantum state discrimination that were 
derived by Holevo and Helstrom ||. 

We now provide simple sufficient conditions on the op- 
timality of the POVM. If the POVM satisfies the ex- 
tremal Eqs. (|^) and (^) and if the following inequalities 
hold: 



A -PiPj > 0, i ^ 1, 
A - ao- > 0, 



(12) 
(13) 



then the POVM is the optimal one that maximizes the 
success rate Ps for a given fixed probability of inconclu- 
sive results Pj. 

To prove this statement we show that the Lagrange 
multipliers provide an upper bound on the success rate 
and that this bound is saturated by the POVM that sat- 
isfies Eqs. (|^) and (||). From the definition of the success 
rate (||), the inequalities ( [l^ ) and the normalization (^) 
we obtain 
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Ps<ETr[^nj]=Tr[A(l-no 



(14) 



Now we use the inequality (nSh and finally we take into 
account the constraint Tr[(Tno] = Pj. We arrive at 



Ps < Tr[A] - aPi. 



(15) 



This last inequality shows that Pg is limited from above 
by the quantity that depends only on the Lagrange mul- 
tipliers A and a and also on the fixed P\. If the POVM 
IIj satisfies the extremal Eqs. (^) and (||) then this upper 
bound is reached, as can easily be checked. 

We have thus established a simple criterion for check- 
ing of the POVM optimality. Of course, we would like to 
derive the optimal POVM 11^ for given pj, pj and P\. The 
analytical solution to this problem seems to be extremely 
complicated. Nevertheless, recently it was pointed out 
that one can solve this kind of problems very efficiently 
numerically p^ . One possible simple and fruitful ap- 
proach is to solve t he extremal e quat ions by means of 



repeated iterations l34 



one could iterate directly Eqs. (Th and 




In principle. 
However, the 



POVM elements 11^- should be positive semidefinite Her- 



mitian operators. All constraints can be exactly satisfied 
at each iteration step if the extremal equations are sym- 
metrized. First we express lij = pjX^^ pjllj and combine 
it with its Hermitian conjugate. We proceed similarly 
also for Ho and we get 



U,=p^\-'p,Il,p,\-\ 
Ho = a^A^VHocrA^^ 



j = l,...,7V, 



(16) 
(17) 
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The Lagrange multipliers A and a must be determined 
self-consistently so that all the constraints will hold. If 
we sum Eqs. (|l^) and ( [l7| ) and take into account that 

EN 
j=0 



Hi = 1, we obtain 



N 



1/2 



(18) 



The fraction of inconclusive results calculated for the 
POVM after the iteration is given by 



Pi = a^Tr[crA"VnocrA"^]. 



(19) 



Since the Lagrange multiplier A is expressed in terms 
of a, Eq. ( [l9|) forms a nonlinear equation for a single 
real parameter a (or, more precisely, a^). This nonlin- 
ear equation can be very efficiently solved by Newton's 
method of halving the interval. At each iteration step for 
the POVM elements, we thus solve the system of coupled 
nonlinear equations (18) and (|l^) for the Lagrange multi- 
pliers. These self-consistent iterations work very well and 
our extensive numerical calculations confirm that they 
typically exhibit an exponentially fast convergence ]34| . 

We note that the maximization of the success rate Pg 
for a fixed fraction of inconclusive results Pi can also be 
formulated as a semidefinite program. Powerful numeri- 
cal methods developed for solving this kind of problems 
may be applied. Here we will not investigate this issue in 
detail and we refer the reader to the papers 

(H i^, 

where the formulation of optimal quantum-state discrim- 
ination as a semidefinite program is described in detail. 
Note also that the semidefinite programming has recently 
found its applications in several branches of quantum in- 
formation theory such as the optimization of completely 
positive maps |4^, Q , the analysis of the distillation 
protocols that preserve the positive partial transposition 
psf , or the tests of separability of quantum states |E6| . 



III. MAXIMAL ACHIEVABLE RELATIVE 
SUCCESS RATE 



As the fraction of inconclusive results is increased the 
success rate Ps decreases. However, the relative success 
rate Prs grows until it achieves its maximum. If {pj} 



N 



are linearly independent pure states, then Pas.max = 1 
because exact IDP scheme works and the unambiguous 
discrimination is possible. Generally, however the max- 
imum is lower than unity. To find this maximum, we 
notice that in the limit Pi — > 1 the POVM element Hq 
must tend to the identity operator. This means that at 
some point Hq becomes positive definite operator and all 
its eigenvalues are strictly positive. In that case, the ex- 
tremal equation (^) can be satisfied if and only if 



A = 



(20) 



Since we are looking for some nontrivial solution to 
the extremal equations with Pi < 1, at least one of the 



extremal Eqs. (0) must have a nontrivial solution Hj ^ 0. 
This implies that at least one of the operators A — PjPj 
must have one eigenvalue p, equal to zero which implies 
that 



det [acr — PjPj] — 



(21) 



must hold at least for one of the states pj. The optimal 
Hj is then proportional to the projector to the subspace 
spanned by eigenvectors corresponding to the eigenvalue 
p = 0. 

The maximal attainable relative success rate is ob- 
tained if wc insert (|2^) into (|l5|), take into account that 
Tr[(T] = ^jPj = 1 and re- normalize according to Eq. (|^), 
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R,S 



(22) 



To determine the maximal Prs we must find the maximal 
a that satisfies Eq. (|2l|). Since a is positive definite it can 
be inverted, and we can equivalently express the maximal 
a as the maximal eigenvalue of a Hermitian matrix. 



flj = j3j max[eig((T ' pjU '")] 
The maximal Prs is equal to the largest aj 



RS.max = maxflj. 
i 



(23) 



(24) 



For qubits, Eq. ( |21| ) leads to quadratic equation for the 
multiplier aj that can be solved analytically, 

{a, - p,f = a)Tr[a^] - 2a,p,Tr[ap,] + p^Tr^^] . (25) 

It turns out that the maximal Prs depends only on 
the a-priori probabilities pj, the purities of the states 
Vj — Tr[/9^] and the overlaps Ojk = Tr[pjPfe]. In this 
context it is worth noting that it was shown recently that 
these parameters of the quantum states can directly be 
measured without the necessity to carry out a complete 
quantum state reconstruction Esl M . 



IV. DISCRIMINATION BETWEEN TWO 
MIXED QUBIT STATES 

We proceed to illustrate the methods developed in 
the present paper on explicit example. We consider the 
simple yet nontrivial problem of optimal discrimination 
between two mixed qubit states pi and p2- To sim- 
plify the discussion, we shall assume that the purities 
of these states as well as the a-priori probabilities are 
equal. Pi = V2 = P, Pi = P2 = 1/2. The mixed states 
can be visualized as points inside the Poincare sphere 
and the purity determines the distance of the point from 
the center of that sphere. Without loss of generality, we 
can assume that both states lie in the xz plane and are 
symmetrically located about the z axis. 



Pi, 2 = vi^iA^) 



1 



1, 



(26) 
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where the parameter 77 determines the purity, 
\ipj){'ipj \ denotes a density matrix of a pure state, 



1^1,2(0)) 



cos-|0)±sin-|l) 



(27) 



and E (0,7r/2). From the symmetry it foUows that the 
elements Hi and 112 of the optimal POVM must be pro- 
portional to the projectors 'ipi{4>) ^md ip2{4')i where the 
angle </) G (7r/2,7r) is related to the fraction of the incon- 
clusive results. The third component XIq is proportional 
to the projector onto state jO). The normalization of the 
POVM elements can be determined from the constraint 
(0) and we find 



ni,2(0) 
no(0) 



1 



2sin^((/)/2) 



1 



^1,2(0), 



tan^(</)/2) 



|0)(0|. 



(28) 



The relative success rate for this POVM reads 

p _ 1 + r]cos{4> - 6) 
2(1 -f 77 cos 6* cos 0) 

and the fraction of inconclusive results is given by 



Pi = i (1 + 77 cos 6*) f 1 }, , , 

2 ^ ' ^ V tan2(0/2) 



(29) 



(30) 



The formulas ( p^ ) and ( pO| ) describe implicitly the depen- 
dence of the relative success rate Prs on the fraction of 
the inconclusive results Pi. From Eqs. (|l^) and (|l^) one 
can determine the Lagrange multipliers A and a for the 
POVM (^and check that the extremal Eqs. (0), (|), 
(12) and (|13|) are satisfied which proves that the POVM 
(28) is indeed optimal. 

The maximum PRs.max (^^ is achieved if the angle 
is chosen as follows, 



C0S( 



On inserting the optimal 



— —77 cos 0. (31) 
max back into Eq. (p9) we get 



77 sm c 



\/\ — rf cos^ 6 



(32) 



Making use of Eqs. (^J) and (^5|) we can express the 
PRS,max in terms of the overlap and the purity V, 
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RS , max 



P-Oi 



2 - p - O12 



(33) 



If we calculate O and V for the density matrices (|2q) a nd 
insert them into ( |3^ ) then we recover the formula p2[). 

The optimal POVM ( p8| ) can be also obtained numer- 
ically. We demonstrate feasibility of iterative solution 
of the symmetrized extremal equations (p^), (p^), (p8[), 
and (^^ for mixed quantum states (|2^) with the angle of 
separation 9 = 7r/4. The trade-off of the relative success 
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FIG. 1: Relative success rate Pas versus the fraction of in- 
conclusive results Pi for the optimal discrimination of two 
mixed states ( p6[ ) with 6 = tt/A and four different parameters 
■q = 0.7, -q = 0.8, 77 = 0.9, and r) = 1.0. 



rate and the probability of inconclusive results is shown 
in Fig. for various purities of the states being discrimi- 
nated. For the given probability Pi of inconclusive results 
and the given purity of the states the extremal equations 
are solved self-consistently by means of repeated itera- 
tions. The success rate Prs is evaluated from the ob- 
tained optimal POVM and re-normalized according to 
Eq. (H). The numerically obtained dependence of Prs on 
Pi is in excellent agreement with the a naly tical depen- 
dence following from formulas ( |29| ) and (|3C|). Typically, 
a sixteen digit precision is reached after several tens of it- 
erations. The trade-off curves shown in Fig. 1 reveal the 
monotonous growth of Prs until the maximal plateau 
(^) is reached. 



V. CONCLUSIONS 

In conclusion, we have considered a generalized dis- 
crimination scheme for mixed quantum states. The 
present scenario interpolates between the Helstrom and 
IDP schemes. We allow for certain fixed fraction of incon- 
clusive results and maximize the success rate. We have 
derived the extremal equations for the optimal POVM 
that describes the discrimination device. The extremal 
equations can be very efficiently solved numerically by 
means of the devised simple iterative algorithm or, alter- 
natively, by using the powerful techniques of semidcfinite 
programming. 

We have showed that the relative success rate Prs 
monotonically grows as the fraction of inconclusive re- 
sults Pi is increased and at certain point it reaches its 
upper bound Prs, max- For pure linearly independent 
states this bound is Prs, max = 1 which corresponds to 
the IDP unambiguous discrimination scheme. For mixed 
states this bound is in general lower than unity and we 
have derived a simple formula for it. 
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The present scheme may be important for quantum 
cryptographic schemes where the receiver and/or eaves- 
dropper want to discriminate nonorthogonal states. Al- 
though these schemes are usuahy based on pure states, 
in rcahstic cases the unavoidable noise and decoherence 
will reduce the purity of these states and one would have 
to deal with mixed states. Various modifications of our 
method can be suggested for such application. For in- 
stance, if the involved states are in some sense asymmet- 
ric, one may impose the condition that the probabilities 
of inconclusive results or successful discrimination of the 
states pj should all be identical, and minimize the error 



rate with this additional constraint. 
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